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ABSTRACT 

Context. We present a study of the Silicon-bearing species detected in a line confusion limited survey towards Orion KL performed 
with the IRAM 30-m telescope. 

Aims. The analysis of the line survey is organized by families of molecules. Our aim is to derive physical and chemical conditions 
for each family taking into account all observed lines from all isotopologues of each species. Due to the large number of transitions 
in different vibrationally excited states covered by our data, which range from 80 to 280 GHz, we can provide reliable source average 
column densities (and therefore, isotopologue abundances and vibrational temperatures) for the detected molecules. In addition, we 
provide a wide study of the physical properties of the source based on the different spectral components found in the emission lines. 
Methods. We have modeled the lines of the detected molecules using a radiative transfer code, which permit us to choose between 
Large Velocity Gradient (LVG) and Local Thermodinamic Equilibrium (LTE) approximations depending on the physical conditions 
of the gas. We have used appropriate collisional rates for the LVG calculations. In order to qualitatively investigate the origin of the 
SiS and SiO emissions in Orion KL we ran a grid of chemical models. 

Results. For the v= 1 state of SiO we have detected the J=2- 1 line and, for the first time in this source, emission in the J=4-3 transition, 
both of them showing strong masering effect. For SiO v=0, we have detected 28 SiO, 29 SiO, and 30 SiO; in addition, we have mapped 
the J = 5-4 SiO line. For SiS, we have detected the main species, 29 SiS, and SiS v=l. Unlikely other species detected in Orion KL 
(IRc2), the emission peak of SiS appears at a velocity of =T5.5 km s _1 ; a study of the 5-4 SiO line around IRc2 shows this feature as 
an extended component that probably arises from the interaction of the outflow with the ambient cloud. We derive a SiO/SiS column 
density ratio of a* 13 in the plateau component, four times lower than the cosmic O/S ratio a* 48. In addition, we provide upper limits to 
the column density of several non-detected Silicon-bearing species. The results of our chemical models show that while it is possible 
to reproduce SiO in the gas phase (as well as on the grains), SiS is a product of surface reactions, most likely involving direct reactions 
of Sulphur with Silicon. 

Key words. Surveys - Stars: formation - ISM: abundances - ISM: clouds - ISM: molecules - Radio lines: ISM 



1. Introduction 

The Orion BN/KL (iBecklin & Neugebaueti Il967t 
iKleinmann & Lowl I1967I) nebula is one of the most stud- 
ied star fo rmation regions in t he Milky Way. At a distance 
of 414 pc dMenten et all 120071) the nebula is embedded in a 
giant molecular cloud harboring practically all phases of the 
interstellar medium, from hot and diluted plasma, to PDRs, 
protostellar cores, molecular outflows, SiO and H2O masering 
regions, high density cores, intermediate and high mass star 
formation, protoplanetary disks, and proplydes (see , e.g. , 
Genzel et allll98rAlGenzel & Stozkl [19891: IWright et all 1 1995: 
Cernicharo et allll990lll994HPlambeck et alll2009l) . 



* Appendix A (density diagnostic), Appendix B (online Figures), and 
Appendix C (online Tables) are only available in electronic form via 
http://www.edpscience.org 

** This work was based on observations carried out with the IRAM 
30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain) 



Together with Sgr B2, Orion BN/KL nebula exhibits a rich 
spectrum (see, e.g., Tercero et al., 2010, hereafter Paper I, 
and references therein) produced by complex organic molecules 
which are formed through reactions on the grain surfaces dur- 
ing the collapse phase followed by evaporation when radiation 
from a newly formed star becomes available. Due to the high 
temperature of the gas, molecular lines are particularly strong 
in Orion, allowing several line surveys of this source over the 
last 20 years. Recently, we have performed a line survey to- 
wards Orion IRc2 source between 80 and 280 GHz (Paper I), 
not limited by sensitivity but only by line confusion. The data 
provide a significant number of transitions for all molecules de- 
tected so far towards this source. Although physical structure 
of the Orion is rather complex, the large number of transitions 
observed for each species allows to model the different cloud 
components and to derive reliable physical parameters. In addi- 
tion, the line survey provides a deep insight on the chemistry of 
the Orion KL region and allows to refine our knowledge of its 
chemical structure by searching for new molecular species and 
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new isotopologues and vibrationally excited states of molecules 
already known to exist in this source. In a first paper we have 
presented the line survey and analyzed the CS-bearing species 
deriving the abundance of CS, OCS, CCS, CCCS, H 2 CS, and 
HCS + (Paper I). In this paper we will analyze the Silicon-bearing 
species SiO and SiS. SiO has been observed with sin gle dishes 
and interferometers (see, e.g., Plambe ck et all |2009 and refer- 
ences therein), while only a few observations are available fo r 
SiS (iDickinson & Rodriguez-Kulr^[TMllzi^[T9M[T99TI) . 

SiO lines in Orion show a complex pattern of thermal and 
maser emission. The masers seem to arise from a small regi on 
around the radio continuum source I (Chur chwell et all 1 1987b . a 
young star with a very high luminosity without infrared co unter- 
part, ^10 5 L , dGezari et all fl998t iGreenhill et all l200l . This 
source is also driving the low velocity outfl ow observed in SiO 
dBeuther et all 2005; Plambe ck et all l2009) and in other molec- 
ular species (Paper I). Recent studies have di scussed the driving 
source of the high velocity outflow: whereas [Beuther & Nissenl 

(2008) claimed that SMA1 (a sub-millimeter so urce not detected 
at IR and centimeter w avelengths, predicted bv lde Vicente et all 
120021 and detected by iBeufher et ail |2004|) is the host of the 
high velocity outflow (based on combined observations of 
7=2-1 C 18 Q from the SMA and the IRAM 30m telescope), 
Plambec k et al.l (120091) defended that this outflow is a continua- 
tion of the low velocity outflow (based on CARMA observations 
ofSiOv=0 7=2-1). 

The SiO v=l maser emission was modeled in early inter- 
ferometric obs ervations as ari sing from a rotating and expand- 
ing disk (Plambeck et all fl990t) . H owever, maser emis sion was 
also found in the v=0 7=2-1 line bv lWright etal.1 d 1995b adding 
more complexity to the modeling of th e structure of the emitting 
source. Recently, Godd i et al.l (12009) have found maser emis- 
sion in the v=0 7=1-0 transition of 29 SiO and 30 SiO associ- 
ated w ith source I. In addition, observations of Plambec k et al.l 

(2009) clearly show that the emission of this line arise essentially 
from an outflow driven by source I. Although our single dish 
data cannot provide a view of the spatial structure of the ther- 
mal and maser emission around that source, the observed maser 
and thermal lines can provide useful constraint on the physical 
conditions of the gas. 

In addition to the study of SiO and SiS we derive upper limits 
to the abundance of SiC, SiC 2 , c-SiC 3 , SiC 4 , SiN, SiCN, SiNC, 
ob-SiC 3 , l-SiC 3 , Si 3 , SiCCO, SiCCS, SiH 2 , H 2 CSi, and the dif- 
ferent isomers of Si 2 H 2 . The observations are described in Sect. 
[2] The results for SiO and SiS are analyzed in Sect. [3] Section 
3]is devoted to the modeling of the observed lines. All these re- 
sults are discussed in Sect.|5]in terms of comparisons with chem- 
ical models predictions for Silicon-bearing species. The effect 
of the new collisional rates for SiO and SiS is analyzed in the 
Appendix. 

2. Observations and Data Analysis 

The observations were carried out using the IRAM 30m radiote- 
lescope during 2004 September (3 mm and 1.3 mm windows), 
2005 March (full 2 mm window), 2005 April (completion of 3 
mm and 1.3 mm windows), 2007 January (maps and different 
positions) and March 2008 (2-D line survey using a multibeam 
receiver). Four SIS receivers operating at 3, 2, and 1.3 mm were 
used simultaneously with image sideband rejections within 20- 
27 dB (3 mm receivers), 12-16 dB (2 mm receivers) and 13 
dB (1.3 mm receivers). System temperatures were 100-350 K 
for the 3 mm receivers, 200-500 K for the 2 mm receivers and 
200-800 K for the 1.3 mm receivers, depending on the particu- 



Table 1. tjmb and HPBW along the covered frequency range 



Frequency (GHz) 


T\MB 


HPBW (") 


86 


0.82 


29.0 


110 


0.79 


22.0 


145 


0.74 


17.0 


170 


0.70 


14.5 


210 


0.62 


12.0 


235 


0.57 


10.5 


260 


0.52 


9.5 


279 


0.48 


9.0 



Notes. r\ MB and HPBW along the covered frequency range. 



lar frequency, weather conditions, and source elevation. For the 
spectra between 172-178 GHz, the system temperature was sig- 
nificantly higher, 1000-4000 K, due to proximity of the atmo- 
spheric water line at 183.3 GHz. The intensity scale was cali- 
brated using two absorbers at different temp eratures and using 
the Atmospheric Tr ansmission Model (ATM. ICernicharol 1 1 98 5 : 
iPardo et all 12001 al) . The half power beam width (HPBW) and 
the main beam efficiency (t\mb) along the covered frequency 
range are given in Table Q] 

Pointing and focus were regularly checked on the nearby 
quasars 0420-014 and 0528+134. The observations were made 
in the balanced wobbler-switching mode, with a wobbling fre- 
quency of 0.5 Hz and a beam throw in azimuth of ±240". No 
contamination from the off position affected our observations 
except for a marginal one at the lowest elevations (~ 25°) for 
molecules having emission along the extended ridge. 

The backends used were two filter banks with 512x1 MHz 
channels and a correlator providing two 512 MHz bandwidths 
and 1.25 MHz resolution. We pointed towards the (survey) posi- 
tion a = 5 A 35'" 14.5 s , 5 = -5° 22' 30.0" (J2000.0) corresponding 
to IRc2. 

The spectra shown in the figures are in units of main beam 
temperature, Tmb- In spite of the good image band rejection of 
the receivers, each setting was repeated at a slightly shifted fre- 
quency (10-20 MHz) in order to identify and remove all features 
arising from the image side band. In the data reduction we have 
removed most of them above 0.05 K (see Paper I for further ex- 
planation of this procedure). 

2.1. 2D survey observations 

The HERA multibeam receiver was used to carry out a sys- 
tematic line survey between 216-250 GHz over a region of 
144x144" (2x2' approximately) with a 4" spacing. System tem- 
peratures were around 400-500 K. The sensitivity of the result- 
ing maps was comparable to those of the pointed line survey 
with the single pixel SIS receivers. The 2D line survey obser- 
vations were performed in On-The-Fly mapping mode, with po- 
sition switching using a emission-free reference position at an 
offset (-600",0) from IRc2. We used WILMA as the main back- 
end, covering the full 1 GHz bandwidth provided by HERA with 
2 MHz of spectral resolution. We also used in parallel the ver- 
satile VESPA spectrometer to get some interesting lines within 
the 1 GHz range with higher spectral resolution (320 kHz), as 
we did for 7=5^1 SiO (from all the 2D survey data only the map 
of this line is showed in this paper). A full description of the 
2D line survey will be published elsewhere (Marcelino et al., in 
preparation). 
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3. Results 

The line survey has been presented elsewhere (Paper I). This pa- 
per is devoted to the study of the emission of Silicon-bearing 
molecules: we have detected SiO, 29 SiO, 30 SiO, SiO v=l, SiS, 
29 SiS, and SiS v=l. For Si ls O and Si 17 we have evidences 
of their presence but all their lines in our frequency range are 
blended with other molecules. We report the detection of the 
transition J=4-3 of SiO v=l for the first time in this source. 
While this transitions shows strong maser emission, rotational 
lines for v>2 in the covered frequencies lack of maser emission. 
In the following subsections we describe the data and the results 
for each molecular species. 

3.1. SiO 

Five rotational transitions of silicon monoxide (SiO) fall in the 
covered frequency range of our line survey. For its isotopologues 
( 29 SiO, 30 SiO, Si ls O, and Si n O) and SiO vibrationally excited 
four rotational transitions are within our frequency range. 

The rotational constants used to derive the line parame- 
ters have been taken from ISanz et all (I2003I) (for SiO, 29 SiO, 
30 SiO, Si ls O, SiO v=l) and for Si n O, the line parameters have 
been derived from the Dunham coefficients, YyS, of 28 SiO and 
the isotopic relations Y,/=Y, 7 *(//(Si 17 0)/Ai(SiO)) (/+2 -'' )/2 (where 
ytz is the reduced mass of the i sotopologue). The dipo le moment 
(/z=3.098D) was reported by Ravmon da et alJ dl970l) . Line pa- 
rameters and observed intensities for all these lines are given in 
Table[2] FigureQ]shows the observed lines of SiO and its isotopo- 
logues, together with the results from our best model (see Sect. 
I4.1l i. Figure [2] shows the SiO v=l lines within our line survey. 
The 7 = 2-1 and 7 = 4-3 lines show maser emission whereas the 
7 = 5-4 and 7 = 6-5 are blended with other abundant molecules 
in Orion and their possible contribution to the observed features 
is rather weak and two orders of magnitude below the observed 
emission in the 7=2-1 and 7=4-3 maser emission. Rotational 
lines from SiO v >2 are below the detection limit of this line 
survey. Only the v=2 7=2-1 seems to have a very week emis- 
sion, Tmb -0.06 K, at the velocity of the red component of the 
maser. 

The line profiles indicate the contribution to the emerging 
intensities from different velocity components: the extended and 
compact molecular ridge (difficult to separate), the hot core, the 
high velocity outflow, and the low velocity plateau. The latter 
appears as the most significant contribution to the line profile 
and intensity. 

Table IC.ll available electronically only, gives the parame- 
ters for the different cloud components derived by fitting several 
Gaussian profiles to the observed transitions of SiO and to se- 
lected lines of 29 SiO and 30 SiO. The profiles are separated into 
two components: one wide component which corresponds to the 
low and high velocity outflows (the plateau), and a narrow one 
consisting of a mix of the ridge components with the plateau. 

The broad component from the plateau dominates the line 
profiles and hides the contribution of the hot core at 5 km s _1 , 
which is only marginally detected as a blue shoulder in the emis- 
sion from the isotopologues. We interpret this behavior as due to 
the opacity of the SiO lines in the high velocity gas which could 
absorb the e mission from the h ot core (this effect was studied 
for HDO by iPardo et al.ll2001bl) . However, the line profiles of 
the 29 SiO and 30 SiO isotopologues, for which the high velocity 
gas should be thinner, clearly need a component at the velocity 
of the hot core (see also the discussion for SiS below) in order to 
reproduce their line profiles. 



3.1.1. 7=5-4 SiO map 

From the 2D survey data of Orion KL, a map of the SiO 7=5-4 
intensity at different velocity ranges is shown in Fig. IB. II The 
velocity structure of the SiO emission shows the contribution 
from all the cloud components quoted above. Note the spatial 
displacement of the emission peak with velocity. Particularly 
interesting is the spatial distribution of the red and blue wings 
at the largest velocities (panels top left and bottom right). The 
high velocity o utflow appears a s an ell iptical shell of gas around 
IRc2. Recently, Plambeck et al. (2009) have obtained a map with 
an angular resolution of ^0.5" of the v=0 7=2-1 line of SiO. 
Their data for the extreme red velocities, 29-52 kms -1 , indicates 
that the emission is shifted towards the East from source / by 
5-10". However, blue extreme velocities, -20 to -4 km s -1 , are 
found several arcseconds W and NW of source /. IWright et al.l 
dl996 ) obtain a similar result in their aperture synthesis (^4" 
resolution) v=0 7=2-1 SiO velocity map of the extreme veloci- 
ties (vlsr<— 11 km s _1 and \isr>29 km s _1 ). These results are 
in very good agreeme nt with our lower angular resolution map. 
iPlambeck et al.l d2009l) found that the bulk of the emission arises 
from a bipolar outflow covering velocities from -13 to 16 km 
s _1 , driven by source / and with an extent of ^6" along the NE- 
SW direction. The large velocity outflow, or extreme velocities, 
seems to be a continuation of the low velocity outflow but less 
spatially structured. This "EW bipo larity" of the S i O hig h ve- 
locity outflow was noted already by Olofs son et alJ ([.198 1). We 
have obtained angular source sizes between 16" for the central 
velocities to 23" for the extreme velocities, assuming emission 
within the half flux level and corrected for the size of the tele- 
scope beam at the observing frequency. The distribution of the 
high velocity gas is very similar to that foun d in the extended 
maser emis sion of H2O at 183.3 GHz found bv lCernicharo et all 
(I1990LI1994 . 

3.2. SiS 

Nine transitions of silicon monosulfide 28 Si 32 S v=0,l and of its 
rare isotopologues are present in the covered frequ ency range. 
Spectroscopic constants are from lSanz et alJ (|2003l). Th e dipol e 
moment (yU=1.730D) has been taken from iHoeft et alJ |[l969). 
Line frequencies and observational line parameters are given in 
Table [3] The observed line profiles and the results from our best 
model (see Sect. I4.2l i are shown in Fig. [3] for the main isotopo- 
logue, and in Fig. lB.2l (for 29 SiS and SiS v=l). The lines from the 
other isotopologues are blended with strong lines arising from 
other molecular species. Only one line (7=9-8 transition) is free 
of blending for Si 34 S. The line intensities arising from the iso- 
topologues 30 SiS and Si 33 S are below the confusion limit of our 
line survey. 

The line profiles of the most abundant isotopologue display 
three components: a wide component that corresponds to the 
plateau, the hot core (clearly seen in the profile of the 7 = 14- 
13 transition), and a narrower one at Vlsr - 15.5 km s -1 . No 
emission has been observed from the ridge component. Line pa- 
rameters are given in Table IC.21 available electronically only (in 
the wide component emission from the plateau and the hot core 
are merged). 

The velocity of the emission peak for the main isotopologue 
is at Vlsr ^15.5 km s _I which coincides with the LSR veloc- 
ity of the red component of the SiO v= 1 maser emission (which 
could be a fortuitous agreement, see Sect. l3.2.TT i. Previous works 
discuss the presence of SiS in Orion KL and its association 
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Table 2. Lines of SiO, SiO isotopologues, and SiO vibrationally excited 



Molecule 


Observed 


Tmb 


7* 

J T MB dv 


Trasition 


Rest 


E„/k 


Sy 




Wlsb (km s" 1 ) 


(K) 


(Kkms- 1 ) 


J 


Frequency (MHz) 


(K) 




SiO 


11.0 


17.5 


394+18 


2-1 


86846.971 


6.3 


2.00 




8.5 


25.6 


590+18 


3-2 


130268.665 


12.5 


3.00 




8.5 


44.6 


962+29 


4-3 


173688.210 


20.8 


4.00 




7.3 


45.0 


1485+22 


5-4 


217104.889 


31.3 


5.00 




8.0 


61.1 


2219+43 


6-5 


260517.985 


43.8 


6.00 


29 SiO 


11.5 


1.27 




2-1 


85759.202 


6.2 


2.00 




11.3 


9.22 


238+20 


4-3 


171512.814 


20.6 


4.00 




15.5 1 


9.40 




5-4 


214385.778 


30.9 


5.00 




10.6 2 


12.9 


... 


6-5 


257255.249 


43.2 


6.00 


30 SiO 


14.9 3 


1.41 




2-1 


84746.193 


6.1 


2.00 




10.1 


6.17 


108+5 


4-3 


169486.929 


20.3 


4.00 




10.6 


6.36 


140+4 


5-4 


211853.546 


30.5 


5.00 




10.0 


12.6 4 


153+14 


6-5 


254216.751 


42.7 


6.00 


Si ls O 


6.9 5 


0.26 




2-1 


80704.907 


5.8 


2.00 




8.6 


4.30 5 - 6 




4-3 


161404.866 


19.4 


4.00 




11.2 


0.65 2 




5-4 


201751.443 


29.0 


5.00 




7 






6-5 


242094.928 


40.7 


6.00 


Si n O 


20.4 8 


0.27 




2-1 


83588.641 


6.0 


2.00 




12.3 


0.40 




4-3 


167171.974 


20.1 


4.00 




12.7 


0.7 1 9 




5-4 


208959.990 


30.1 


5.00 




9.2 


0.55 




6-5 


250744.688 


42.1 


6.00 


SiO v = 1 


14.2 
-3.9 


77.5 
113.5 


587+20 
750+19 


2-1 


86243.430 


6.2 


2.00 




14.0 


21.5 


123+14 


4-3 


172481.125 


20.7 


4.00 




-4.5 


4.20 














10 


<1.2 




5-4 


215596.029 


31.0 


5.00 




2 


<1.2 




6-5 


258707.350 


43.5 


6.00 



Notes. Emission lines of SiO, SiO isotopologues, and SiO vibrationally excited present in the frequency range of the Orion KL survey. Column 1 
indicates the species, Col. 2 gives the observed (centroid) radial velocities, Col. 3 gives the peak line temperature, Col. 4 the integrated intensity. 
Col. 5 the line transition, Col. 6 the calculated rest frequencies, Col. 7 the energy of the upper level, and Col. 8 gives the line strength. 
(1) Blended with 34 SH 2 . (2) Blended with CH 3 OCOH. (3) Blended with CH 3 OH. <4) Blended with CH 3 CH 2 CN b type. (5) Blended with CH 3 CH 2 CN 
in the plane torsion (Tercero et al., in preparation). <6) Blended with CH 2 CHCN. <7) Blended with CH 3 CH 2 CN a type. (8) Blended with H 53 /3 and U 
line. (9) Blended with CH 3 CH 2 C 15 N. (10) Blended with CH 3 CH 2 CN out of plane torsion (Tercero et al., in preparation). 



with the SiO v-1 maser ( Dickinson & Rodrfguez-Kuiperl 198 lb 



ISutton et all fT^l IZiurvsL 1 19881 11991b LSchilke et al., 1997). 



Although confusion is large when considering the weak lines 
of 29 SiS and SiS v=l, we found that their emission peaks at Vlsr 
=* 13.5 km s (see Fig. IB. 21 ). a mixture of all cloud components, 
but dominated by the 15.5 km s feature. 



3.2.1 . The feature at 1 5.5 km s" 1 

At the position of the line survey (IRc2) only SiS (v=0, 1) lines 
and one component of the Si O maser emission ( v=l) show an 
intensity peak at 15.5 km s -1 . IWright et alj (jl990) measured the 
absolute position of the SiO masers to an accuracy of 0".15 co- 
inciding with the position of radio source / (very close to IRc2). 
In order to check the origin of this feature we have observed the 
line profiles of several molecules around IRc2. We have found 
that the line profiles of SiO show important changes with posi- 
tion. Figure [4] central panel displays a 20" x 20" map centered 
on IRc2 of the 7=5-4 line of SiO; the small panels around the 
map show the line profile of the 7=5-4 transition of SiO at se- 
lected positions and show strong differences across the cloud. 
The intensity peak of SiO lines is at 15.5 km s at the positions 
Aa=-10", A6-+6" and Aa=-6", A<5=+14". A comparison of the 
7=5-4 SiO line at (-10", 6") with the emission of the 7=14-13 
SiS line at (0, 0) is shown in the right panel. 



In our map, this extended feature is seen around the posi- 
tion Aa=-7", A<S=+7" (near the BN object) with a radius of 
^5 arcseconds. Figure IB. 31 shows emission lines from differ- 
ent molecules at a position offset (-15", 15") from IRc2. Only 
molecules with a strong emission from the plateau component 
(SO, 34 SO, S0 2 , and SiO) show the 15.5 km s-1 component 
at this position. Hence, this feature is not a particularity of SiS 
but arises probably from the interaction of the outflow with the 
ambient cloud. Consequently, lines toward Orion KL showing a 
strong intensity and eventually high opacity in the plateau could 
hide interesting details of other components (hot core, feature at 
15.5 km s-1), making line interpretation a very difficult task. 

4. Physical parameters of the clouds 

Column densities for all detected species have been calculated 
using an excitation and radiative transfer code developed by J. 
Cernicharo (Cernicharo 2010, in preparation). Depending on the 
selected molecule or ph ysical conditions, we assume the large 
velocity gradient (LVG: ISobolevl 1 1 958b ISobolevL 1 1 960h or local 
thermodynamic equilibrium (LTE) approximations. Table [4] re- 
sume the physical parameters we have obtained for each spectral 
cloud component from the modeling of SiO and SiS (note that 
for the new feature at 15.5 km s we have derived the follow- 
ing parameters from the modeling of the SiS lines: T^=200 K, 
n(H 2 )=5x!0 6 cm 4 , and v /m// power ,> !fe .„. 5 ,Y V =7.5 km s" 1 ). We as- 
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J=2-i 



J=4-3 



J=5-4 



J = 6-5 



Tub (K)F^S 






"SiO 




Si 




-100 -50 50 100 
V L3H ( km s_1 ) 



1 50 
V LSE ( km s ~ ) 




V LSR (km s ) 



00 -50 50 100 
V LSR ( km s_1 ) 



Fig. 1. Observed lines (histogram spectra) and model (thin curves) of SiO and its isotopologues. The dashed line shows a radial 
velocity at 9 km s _1 . 



sume uniform physical conditions: kinetic temperature, density, 
radial velocity, and line width. We adopted these values from the 
data analysis (Gaussian fits and an attempt to simulate the line 
widths and intensities with LTE and LVG codes) as representa- 
tive parameters for the different species. Our modeling also takes 
into account the size of each component and its offset position 
with respect to IRc2. Corrections for beam dilution are applied 
for each line depending on their frequency. 

The only free parameter is, therefore, the column density of 
the corresponding observed species. Taking into account the re- 
duced size of most cloud components the contribution from the 
error beam is negligible except for the extended ridge which has 
a small contribution for all observed lines. 

In addition to line opacity effects, we discussed other sources 
of uncertainty in Paper I. 

4.1. SiO 

LTE conditions have been assumed for the hot core, while LVG 
calculations have been performed for the extended ridge, com- 
pact ridge, plateau, and high velocity plate au, with collisional 
cross sections of SiO-p-H2 taken from lDavou & B alanca (2006) 
(see Appendix) for temperatures between 10 K to 300 K includ- 
ing levels up to J = 20 (E M/) =437 K). 

SiO lines appear to be optically thick and therefore the de- 
rived SiO column density could be significantly underestimated. 
The lines of the isotopologues are, however, mostly optically 
thin so we can estimate the SiO column density assuming stan- 
dard isotopic abundance ratios ( 28 Si/ 29 Si-20 and 28 Si/ 30 Si-30; 



And ers & Grevessd fl989). The results obtained are shown in 
Table|5] Due to the weakness of the less abundant isotopologues 
lines (Si 18 and Si I7 0) and large line overlap problems, we can 
only obtain upper limits for their column density. We have esti- 
mated the uncertainty between 20-30 % for the 29 SiO and 30 SiO 
results. 

The cloud component with the largest column density corre- 
sponds to the plateau with N(SiO) p i ateau = (4.7±1.0)xl0 15 cirT 2 ; 
while the total column density of SiO is JV(SiO)=(7.4+2.0)xl0 15 
cm" 2 , both results in agreement with those o btained from the lin e 
survey of Orion using the Odin satellite dPersson et a n i2007h . 
Our study provides a column density between six times to one 
order of magnitu de larger than that obtain ed in the surveys a t 
high frequency by Schil ke et al.l d2001l) and lComito et alJd2 005). 
Larger colu mn densities have be en repo rted bvlZiurvs & Friberg 
( 19871) and IWright et alJ (1 19961) . while iJohansson et all (11984 
and Sutto n et al.ldl995l) obtained a beam average column density 
in good agreement with our results. These differences are mostly 
due to the different assumptions on the physical conditions and 
cloud structure made in the interpretation of the observations. 

We derive a column density 29 SiO/ 30 SiO ratio of 2, 2, 1.7, 
1.4, 1.25 for the extended ridge, compact ridge, plateau, high 
velocity plateau, and hot core, respectively , in agreement with 
the st andard solar value of 29 Si/ 30 Si ^1.5 dAnders & Grevessa. 
1989). We can estimate a lower limit for the isotopic ra- 
tio 16 0/ 18 and an estimation of 18 0/ 17 by means of the 
total column density ratios N(Si 16 0)/W(Si 18 0) £ 239 and 
7V(Si 18 0)/7V(Si 17 0) - 2 with large uncertainty due to the se- 
vere blending of their lines with other species. These values are 
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Table 3. Lines of SiS, 29 SiS, and SiS v=l 



Molecule 


Observed 


Tmb 


J T MB dv 


Trasition 


Rest 


E„/k 


s« 




Vlsr (kms* 1 ) 


(K) 


(Kkms _1 ) 


J 


Frequency (MHz) 


(K) 




SiS 


12.5 


0.29 


6.1+0.3 


5-4 


90771.558 


13.1 


5.00 




16.7 


0.53 


11.5+0.5 


6-5 


108924.294 


18.3 


6.00 




15.5 


1.34 


30±1 


8-7 


145227.045 


31.4 


8.00 




16.5 


2.14 1 


43±1 


9-8 


163376.773 


39.2 


9.00 




16.2 


2.42 




11-10 


199672.219 


57.5 


11.00 




11.1 


3.10 




12-11 


217817.651 


68.0 


12.00 




8.9 2 


20.0 




13-12 


235961.363 


79.3 


13.00 




17.7 


4.45 




14-13 


254103.213 


91.5 


14.00 




15.2 


4.25 




15-14 


272243.058 


104.5 


15.00 


29 SiS 


16.6 


0.019 




5-4 


89103.781 


12.8 


5.00 










6-5 


106923.019 


18.0 


6.00 




16.1 


0.087 




8-7 


142558.873 


30.8 


8.00 










9-8 


160375.213 


38.5 


9.00 
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12-11 


213816.228 


66.7 


12.00 




13.1 


0.25 




13-12 


231626.772 


77.8 


13.00 
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14-13 


249435.521 


89.8 


14.00 




12.9 


0.24 




15-14 


267242.338 


102.6 


15.00 


SiS v=l 


noise level 






5-4 


90329.891 


13.0 


5.00 




2 






6-5 


108394.292 


18.2 


6.00 




12.4 


0.12 




8-7 


144520.370 


31.2 


8.00 




13.3 


0.46 7 




9-8 


162581.760 


39.0 


9.00 




11.9 


0.17 




11-10 


198700.526 


57.2 


11.00 
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12-11 


216757.615 


67.6 


12.00 




14.8 


0.31 




13-12 


234812.983 


78.9 


13.00 




13.8 


0.36 




14-13 


252866.487 


91.0 


14.00 




i 






15-14 


270917.984 


104.0 


15.00 



Notes. Emission lines of SiS, 2, SiS, and SiS v=l present in the frequency range of the Orion KL survey. Column 1 gives the species, Col. 2 the 
observed (centroid) radial velocities, Col. 3 gives the peak line temperature, Col. 4 the integrated intensity, Col. 5 the line transition, Col. 6 the 
calculated rest frequencies, Col. 7 the energy of the upper level, and Col. 8 gives the line strength. 

<" Blended with CH 3 OCOH. < 2 > Blended with 13 CH 3 OH. < 3 > Blended with U line. < 4 > Blended with (CH 3 ) 2 CO. (5) Blended with 34 S0 2 . (6) Blended 
with CH 3 OH. (7) Blended with 33 S0 2 . <8) Blended with CH 3 CH 2 CN b type. 

Table 4. Physical parameters of Orion KL components 



Parameter Extended ridge Compact ridge Plateau High velocity plateau Hot core 15.5 km s 1 component 



Source diameter (") 


120 


15 


30 


30 


10 


10 


Offset (IRc2) (") 





7 








2 


5 


n(H 2 ) (cm- 3 ) 


l.OxlO 5 


l.OxlO 6 


1.0x10 s 


l.OxlO 6 


5.0xl0 7 


5.0xl0 6 


TV (K) 


60 


110 


125 


125 


225 


200 


Vhalf power intensity (km S ) 


4 


4 


25 


50 


10 


7.5 


vlsr (kms' 1 ) 


9 


7.5 


9 


9 


5.5 


15.5 



Notes. Obtained physical parameters for Orion KL. 
Table 5. Column densities - SiO 



Species 


Extended ridge 


Compact ridge 


Plateau 


High velocity plateau 


Hot core 




NxlO 15 (cm" 2 ) 


N xlO 15 (cm- 2 ) 


NxlO 15 (cm- 2 ) 


N xlO 15 (cm- 2 ) 


N xlO 15 (cm- 2 ) 


SiO from 29 SiO 


0.04±0.01 


0.20±0.05 


5.0±1.0 


1.4+0.4 


1.0+0.3 


SiO from 30 SiO 


0.03±0.007 


0.15±0.04 


4.5±1.0 


1.5+0.4 


1.2+0.3 


SiO (Average) 


0.035±0.009 


0.17±0.05 


4.7±1.0 


1.4+0.4 


1.1+0.3 


29 SiO 


0.0020±0.0005 


0.010±0.003 


0.25±0.06 


0.07±0.02 


0.05±0.02 


30 SiO 


0.0010±0.0003 


0.005+0.002 


0.15±0.04 


0.05±0.02 


0.04±0.01 


Si 18 


<0.0005 


<0.0005 


<0.01 


<0.01 


<0.01 


Si 17 


<0.0001 


<0.0005 


<0.005 


<0.005 


<0.005 



Notes. Column densities of the different SiO isotopologues calculated with LVG and LTE codes. 
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-50 v LSR (km s" 1 ) 50 



Fig. 2. Observed lines of vibrationally excited SiO v = 1 showing maser emission in the transitions 7 = 2-1 and J = 4-3. The dashed 
lines show radial velocities at at -3.4, +5.1 and +13.9 km s . 



two times lower than those of the solar system. In Paper I we 
derive d 7V( 16 OCS)//V( 18 OCS) = 250+135, while iPerssonetal] 
(120071) obtained a 18 0/ 17 ratio of 3.6 from C 18 0/C 17 0. 

We have also modeled the 5-4 SiO line in the offset posi- 
tion (-10",6") from IRc2 (see Fig. |4j> in order to provide a col- 
umn density for the 15.5 km s _1 feature. Using the same col- 
umn densities obtained above for the different components, we 
have added the feature at 15.5 km s _1 with a column density of 
(1.0±0.3)xl0 15 cirT 2 . 

4.1 .1 . SiO maser lines 



SiO m aser emission in Orion was discov ered by Snyder & Buhll 
(1 1974 through its v=l 7=2-1 transition. IPlambeck et al.1 (1 19901) 
suggested that the emission arises from an expanding rotating 
disk around IRc2. The emission from this line, studied by 
several authors, has been found to be confined to a region 



of 20-10 AU (3-10xl0 14 cm) around the radiocontinuum 
source / (IChurchwelletall fl987l) and to trace a protostellar 
wind and/or ou t flow expanding with a veloc ity ^ 20 kms 



( Plambec k et all Il990t iMenten & Reidl 119951: IGreenhill et al 



1998 : iDoeleman et all 119991: IPlambeck et all 1 2003 



Doeleman et all l2004t IGreenhill et all l2004r TMatthews et al 



2007 ). Source / has not an infrared counterpart in the range 4-22 
/urn (IGreenhill et ail 120041) . suggesting that it is surrounded by 
a dusty disk producing severe extinction. Source / is placed 
between Orion IRc2 and the hot core and excites the masers o f 
SiO and H z O dGenzel & StuzkiL [19891: IPlambeck et all [1991 . 
Several other lines of SiO in v=l and 2 have been searched 
towards this source an some of them have been detected (v=l 
7=1-0, 7=2-1, 7=3-2; see references above). The v=l 7= 4-3 
line was searched without success bv lSchwartzet al.1 (Il982h . In 
our line survey we have covered the 7=2-1, 4-3, 5-4, and 6-5 
lines of all vibrational states of SiO. As shown in Fig. [2] only 
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the v=l J -2-1 and 4-3 transitions have been unambiguously 
detected (the 7=4-3 for the first time). For the rotational lines 
of the v=2 level the range of velocities covered by the SiO 
v=l J =2-1 maser is always contaminated by other lines. The 
v=2 7=2-1 line shows a weak feature coincident with the red 
emission at 15 km s but it arises from the 42,3-3 12 line of 
E-CH3OCOH as confirmed from the observed intensities of 
other lines of this species. The v=2 7=4-3 line shows a strong 
feature at the velocity of the blue component of the maser but it 
arises from (CH3)2CO. The other lines from v=2, 3, and 4 are 
always blended with other lines. In our survey we have covered 
the same rotational and vibrational quantum numbers of the SiO 
isotopologues. None of them shows any maser effect within the 
sensitivity of our observations. 

The observation of several ro-vibrational lines of SiO can be 
used to constraint the physical properties of the gas at spatial 
scales from a few tens to a few hundred AU from the source. 
iGoddi et all (§009) have modeled the emission of several iso- 
topologues of SiO. They conclude that while the v=0 7=1-0 
maser emission arises from a region with radius >100 AU and 
densities <10 7 cm 3 , the v=l ,2 7=1-0 and 7=2-1 are produced 
in a region of ^100 AU around source /. Our observed inten- 
sity ratio v=l 7=2-1/7=4-3 is ^5-10 depending on the veloc- 
ity component. We obtain such a behavior for column densities 
of ^10 19 cm -2 and densities = 10 8 cm 4 . For higher densities 
the ratio increases and the predicted intensities of the 7=2-1 and 
7=4-3 lines, although large, are not compatible with the inferred 
brightness temperatures which could be well above 10 5 K. We 
have run several models with different physical parameters. Our 
results suggest that the gas responsible for maser emission could 
have densities up to 10 8 cm 4 , sizes of 100 AU, and g as tempera- 
tures around 600-800 K, a factor two below those of Godd i et aD 
(2009). With these parameters for the emitting region the abun- 
da nce of SiO should be as high as 10~ 4 . 

iGonzalez- Alfonso & Cernicha ro ( 1997) have shown that the 
masers of SiO and their isotopologues are affected by line over- 
lap in the infrared between themselves. Such overlaps seem to 
have little effects on the emission of SiO, 29 SiO, and 30 SiO in 
v^O for which maser em i ssion has been found in evolve d stars 
dCernicharo et al.Lll991btlGon zalez-Al fonso et al.Lll996l) . 



4.2. SiS 

In the same way as for SiO, LTE conditions have been assumed 
for the hot core, while LVG calculations have been performed for 
the 15.5 km s 4 feature and the plateau, with collisional cross 
sections S1S-0-H2 fo r 41 levels and 5<Tt<300 K taken from 
iKlos & Liqud (120081) (see Appendix). The fits are shown in Figs. 
[3land lBT2l The column density results are shown in Table[6] We 
have estimated the uncertainty to be about 20-30 % for the re- 
sults of SiS and 50 % for 29 SiS and SiS v=l. Our largest value of 
the column density for SiS corresponds to the feature at 15.5 km 
s 4 obtaining JV(SiS) 15 . 5 kms -i feature = (7.0+1.7) x 10 14 cm 4 ; 
the total column density for SiS is N(SiS), n , „i = (1. 3 5+0.4 0) x 
10 15 cm' 2 . A similar v alue was obtained by Ziurvsl ( 1988b and 
IZiurvsl (1199 11) whereas Dick inson & R odrfguez-Kui peil dl98~fl) 
found N(SiS) = (1-2) x 10 13 cm 4 . 

Assuming that SiS emission is optically thin, as indi- 
cated by our calculations, we derive an isotopic abundance of 
28 Si/ 29 Si-14, 35, and 30 for the 15.5 feature, the plateau, and 
the hot core, respectively i. e. close to the solar system value. 
We provide an average value of 28 Si/ 29 Si = 26+10. 

The N(SiO)/jV(SiS) column density ratio ob- 
served in the plateau is ^13, in good agreement with 



iDickinson & Rodrfguez-Kuipeil (Il98l|) and four time s lowe r 
than the cosmic O/S ratio of 48 (lAnders & GrevesseL Il989h . 
IZiurvsl (11991b found a JV(SiO)/jV(SiS) ratio 40-80. We can 
also derive this ratio by means of Af( 29 SiO)//V( 29 SiS), obtaining 
a value ^ 25 for the plateau. In order to compare column 
density ratios, we have to assume that the region of the line 
formation is the same for each molecule and the excitation 
temperature is similar for both species, for that reason we only 
provide this ratio for the plateau component. In the previous 
work of this line survey (Paper I), we derived O/S ratios using 
different species/families of molecules: A^(HCO + )/A^(HCS + ) 
13, W( H 7 ,CO)/Ar(H?,CS) ^ 1 2, and JV(CO)/JV(CS) * 370; in 
addition, iPersson et all (120071) found Ar(H 2 0)/jV(H 2 S) 20 
and jV(H 2 CO)/jV(H 2 CS) a 15. All the obtained values, but 
N(CO)/N(CS) =i 370, show a similar O/S ratio, indicating that 
the different formation paths of different molecules maintain a 
constant ratio O/S in the same particular region of the cloud. 

^From the column density obtained for SiS in the ground and 
the vibrationally excited states, it is easy to estimate a vibrational 
temperature by means of: 



exp 



(-1?) 



fr 



N(SiS v = x) 
N(SiS) 



(1) 



where E Vj is the excitation energy of the vibrational state (Ey=i = 
1077 K, T v u, is the vibrational temperature, f v is the vibrational 
partition function, N(SiS v = x) is the column density of the 
vibrational state and N(SiS) is the column density of SiS in the 
ground state. The vibrational partition function can be approxi- 
mated by 



/y=l +«/>(-§■) + 

+2exp(-^)+exp(-^) 



(2) 



which, for low T,,,;, leads to f v =s 1 . 

We obtain T V! -j = 500+200 K in the 15.5 km s 4 component. 
This value is higher than the kinetic temperature we have as- 
sumed for that component (200 K). This result could indicate 
an inner and hotter emitting region for vibrationally excited SiS, 
suggesting that the excitation temperature varies across the fea- 
ture at 15.5 km s 4 . In Paper I we calculated the vibrational tem- 
perature in the hot core component for OCS v% = 1 and V3 = 1 
obtaining ^210 K and an upper limit of 300 K for CS v = 1. As 
we have indicated in our previous paper, all these results could 
point to radiative pumping effects in the populations of the vibra- 
tionally excited states of these molecules (taking into account the 
magnitude of the collisional rates of these species). 

Higher angular resolution observations are necessary to re- 
solve any possible excitation gradient and temperature profile in 
the feature at 15.5 km s 4 and in the hot core component. 

4.3. Other Silicon-bearing Molecules 

We provide upper limits for the column density of several Silicon 
molecules not detected in our line survey. We have assumed the 
four typical spectroscopic components of Orion KL (hot core, 
extended ridge, compact ridge, and plateau) plus the 15.5 fea- 
ture, and a LTE approximation for all these molecules. Table [7] 
shows the results obtained, the dipole moment of each species, 
and refer ences for the spectroscopic constants. 

SiC.-|Mi!!ai| {1980) predicted SiC to be an abundant form of 
Silicon in dense clouds on basis of gas phase chemistry models. 
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Table 6. Column densities - SiS 



Species 


15.5 km s 1 Component 
N xlO 14 (cm" 2 ) 


Plateau 

N xlO 14 (cm" 2 ) 


Hot core 

N xlO 14 (cm" 2 ) 


SiS 


7.0±1.7 


3.5±0.8 


3.0+0.7 


29 SiS 


0.50±0.25 


0.10+0.05 


0.10±0.05 


SiS v=l 


0.80+0.40 


0.20+0.10 


0.20±0.10 



Notes. Column densities of different SiS species calculated with LVG and LTE codes. 
Table 7. Column density upper limits 



Molecule 


Column density 
SNxlO 14 (cm" 2 ) 


Dipole 
moment (D) 


References 

spectroscopic constants 


SiC 


1.3 


1.600 1 


(2) 


SiC 2 


0.35 


2.393 3 


(4) 


c-SiC 3 


0.13 


4.200 5 


(6) 


SiC 4 


0.04 


6.420 7 


(8) 


SiN 


0.61 


2.560 9 


(10) 


SiCN 


0.31 


2.900" 


(12) 


SiNC 


0.31 


2.000" 


(12) 


ob-SiC 3 


0.16 


2.200 5 


(13) 


l-SiC3 


0.04 


4.800 14 


(14) 


Si 3 


4 


0.350' 5 


(16) 


SiCCO 


0.40 


1.937 17 


(18) 


SiCCS 


0.65 


1.488" 


(19) 


o-SiH 2 


13 


0.075 20 


(21) 


o-H 2 CSi 


1.7 


0.300 22 


(22) 


p-H 2 CSi 


1.4 


0.300 22 


(22) 


mb-Si 2 H 2 


0.34 


^=0.962/^ =0.039 23 


(24) 


o-db-Si 2 H 2 


2.5 


K^O.480 25 


(25) 



No tes. Upper limits for the c o lumn densi ty of non-detected Silico n-be aring molecules in O rion KL. 

(1) lLanghoff & BauschlicheJ dl990h: (2) ICernicharo et alj dl989l) an d Boaev et al. ( 1990); (3) Suen ram et aD dl989h: ( 4) iGottlieb et alJ dl989l) 
andlSuenram et alj d 1989^(5) I Alberts et alj dl990h: (6)lApponi et all dl999bl): (7)lGordon et alj d2000h: (8)lQhishi et alj dl989T) andlGordon et al l 
(2000 ): (9) | Kerkines & MavridisT ffi)05h: (lO)ISaito et al.) dl983l) andlBizzocchi et all d2006l) : (1 lUApDoni et alj d2000h : ( 12)lApDoni et al] < 2000) 
andlMcCarthy et al] d2001l); (13 ) IMcCarthv et al l dl999l) : (1 4) iMcCarthv et al] d2000l) ; (15 ) IVasiliey et al] dl997l): (16) IMcCarthv & Thaddeusl 
d2003ah: (1 7 )E tschwina (2005); (18) Sanz et al. (2005); (19) Botsc hwina et alj d2002h; (20)lGab riel ( 1993); ( 21)lHirota & Ishikawa ( 1999); (22) 
Izuha et alj ( Tl996h : (23) lGrev & Schaefer IIlldl992h : (24) ICordonnier et al.1 dl992l) and ilVIcCarthv & Thaddeusl d2003bl) : (25) lBogev et al.l d 19941) . 



This molecule has been detected bv lCernicharo et all (Il989h to- 
wards the envelope of the red giant star IRC+101216, but not in 
Orion (see ISchilke et al.L 1 19971 for the attempts of detection in 
this source). We have obtained an upper limit for its total col- 
umn density of 1.3 xlO 14 cirT 2 providing an abundance ratio of 
jV(SiO)//V(SiC);>57. 

SiC -i-- We have n ot detected silicon dicarbide in our line 
survey. iTurned dl99ll) reported the detection of this molecule 
in the hot core of Orion KL. However, this result is quite un- 
certain due to the few observed transitions (three) and their 
weakness. We have searched for other lines of S1C2 in our 
line survey and we conclude that the intensities of the lines 
from this molecule are below the confusion limit as many of 
these lines are blended or missing. We obtain an upper limit 
to its column density of 3.5 xlO 13 cm -2 and an abundance ra- 
tio A^(SiO)/A^(SiC2)S21 1. This molecule and its isotopologues 
have a lso been detected towards I RC + 1 1 2 1 6 by Tha ddeus et al.l 
(11984b (SiC ). ICernicharo etalJd 19861) ( 29 SiC and 30 SiC 2 ) and 
ICernicharo et al.l (Il991al) ( Si^CC ). 

c-SiC- lApponi et al] d!999al) detected rhomboidal SiC 3 in 
the expanding envelope of the evolved carbon star IRC+10216. 
For this molecule we have calculated an upper limit to its 
column density in Orion KL of 1.3xl0 13 cirr 2 , and a ratio 
7V(SiO)/7V(SiC 3 );>570. 



SiC 4.- S1C4 was first detected in the space by lOhishi et al .1 
(1 19891) in the envelope of IRC+10216. For this molecule we have 
calculated an upper limit to its column density in Orion KL of 
3.5xl0 12 cm" 2 , and a ratio 7V(SiO)//V(SiC 4 )£21 14. 

SiN- Sil i con n itride has been detected in IRC+10216 by 
Turn er et all (1 1992]) an d in the galactic center cloud SgrB2(M) 
bv lSchilke et al.l (120031) . We have not found SiN above the con- 
fusion limit in Orion; note, however, that some U-lines in our 
survey could be assigned to SiN (see Fig. IB. 41 panels 1, 6 and 9) 
but we consider that the evidences for its presence in Orion are 
not strong enough within the coverage of our survey. We provide 
an upper limit of 6.1xl0 13 cirr 2 , and JV(SiO)/7V(SiN)>121. 

SiCN.- Cyanosilyli dyne was identified in spectra recorded 
toward IRC+10216 bv lGuelin etail J2000!). The calculated up- 
per limit to its column density in Orion KL is 3.1xl0 13 cm -2 , 
deriving the ratio Af(SiO)/7V(SiCN)>240. 

SiNC- The isocyanosilylidyne isomer has a thermody- 
namic stability very similar to that of SiCN, but a slightly 
smaller dipo l e mom ent. It was detected toward IRC+10216 by 
iGuelin et all (120041) . We provide an upper Iimitof3.1xl0 13 cirr 2 
to its column density, and a ratio 7v'(SiO)/A^(SiNC)^240. 

SiH, SiH^.- In our frequency range the re are not transitions 
of SiH, a molecule tentatively detected by Schi lke et al.l (12001 1) 
towards Orion KL, so we cannot assess that detection. The other 
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Silicon molecule detected in the space (in IRC+10216) is SiH4 
(iGolhaber & Betzi 11984 . but this molecule can be only ob- 
served at IR wavelengths. 

ob-SiCi, l-SiC3, S13, SiCCO, SiCCS, 0-S1H2, 0-H2CS1, p- 
H2CS1, mb-Si2H2, o-db-Si2H2-- These molecules have not been 
detected in the space yet. Upper limits for their column density 
are shown in Table [7] 



5. Discussion 

SiO is a key tracer of shocked emission. Many interferomet- 
ric observations show that thermal and maser SiO emission de- 
picte d the low velocity outflow centered in source / in Orion 



KL dBlakeetal.Lll996t I Wright et all 1 19961: iBeuther et al.l| 2005 ; 
iPlambeck et alTl2009b iGoddi et all l2009TlZapata et all 12009). 
In addition, Si O traces many other molecu l ar outflows in dif- 
ferent sources dlimenez-Serra et all l2004t iGibb et aTL [2 007: 
iDe Buizer et"afl l2009t IZapata et all l2009l) jMookerie a et all 



d2007l) have not found emission from SiO at the position of the 
hot molecular core in G34.26+0. 15. This hot core does not have 
a central source but rather it is externally heated, similar to the 
Orion compact ridge, by sho cks, ionization fronts and stellar 
winds. Mooke riea et all ( 20071) pointed out that the lack of SiO 
in this hot core rules out any significant role played by shocks 
in determining the hot core chemistry. Observations of SiO in 
the L1448-mm outflow permit to distinguis h between the shock 
precu rsor and the postshock components ( Jime nez-Serra et ail, 
120051) : they observed an enhancement in the abundances of SiO 
(and another shock tracers) by one order of magnitude in the 
shock precursor component and three orders of magnitude in the 
postshock gas (leading to the broadening of the line profiles) , 
evide nce of recent ejection of SiO from grains dFlower et all 
1996). 



5. 1 . Molecular abundances 

Molecular abundances were derived using the H2 column den- 
sity calculated by means of the C ls O column density (1.5xl0 16 , 
1.5xl0 16 , lxlO 17 , 5xl0 16 cm" 2 , and 2xl0 17 cirT 2 for the ex- 
tended ridge, compact ridge, plateau, high velocity plateau, and 
hot core, respectively) and the isotopic abundance 16 O/ 18 O=250, 
both provided in Paper I, assuming that CO is a robust tracer 
of H2 and therefore their abundance ratio is roughly constant, 
ranging from CO/H2 - 5xl0" 5 (for the ridge components) to 
2xl0" 4 (for the hot core and the plateau). In spite of the large 
uncertainty in this calculation, we include it as a more intu- 
itive result for the molecules described in the paper. We ob- 
tained N(H 2 ) = 7.5xl0 22 , 7.5xl0 22 , 2.1xl0 23 , 6.2xl0 22 , and 
4.2xl0 23 cm" 2 for the extended ridge, compact ridge, plateau, 
high velocity plateau, and hot core, respectively; for the 15.5 km 
s" 1 component we assume Af(H2)=1.0xl0 23 crrT 2 as an average 
value in Orion KL. In addition, we assume that the H2 column 
density spatially coincides with the emission from the species 
considered. Our estimated source average abundances for each 
Orion KL component are summarized in Table IC.3I (only avail- 
able online), together with comparison valu es f rom other au - 
thors dSutton et al.lll995llPersson et al.ll2007l and|Zjuryj[l988). 
The differences between the abundances shown in Table IC.31 are 
mostly due to the different H2 column density considered, to the 
assumed cloud component of the molecular emission and dis- 
crepancies in the sizes of these components. 



5.2. On the origin of the SiO and SiS emission 

In order to qualitatively investigate the origin of the SiS and SiO 
emissions in Orion KL we ran a grid of m od els using the chemi - 
cal model UCL.CHEM dViti et alll2004allb1:lLerate et alll2010l) . 
a time and depth dependent gas-grain model. 

We modeled the hot core and the plateau separately. Both 
models are two phase calculation. In Phase I we follow the 
chemical and dynamical evolution of a collapsing core up to a 
final density of 5xl0 7 crrT 3 for the hot core component and 10 6 
cm" 3 for the plateau component as derived in Sect. [4] The ini- 
tial gas is at typical densities of ~ 200 cm" 3 and in atomic form 
(apart from a fraction of hydrogen which is alr eady in molecu- 
lar fo rm) and it undergoes a free-fall collapse (Rawlings et all 
1992) until the final densities are reached. During this time, 
atoms and molecules from the gas freeze on to the dust grains 
and they hydrogenate where possible. Note that the advantage 
of this approach is that the ice composition is not assumed but 
it is derived by a time dependent computation of the chemical 
evolution of the gas/dust interaction process. However, the ice 
composition does depend on the percentage of gas depleted on 
to the grains during the collapse, and this in turns depend on the 
density as well as on the sticking c oefficient and other pro perties 
of the species and of the grains (see Rawlings et al .119921) . In our 
model we can vary such percentage (reflecting the uncertainty 
on the grain properties and sticking probabilities) and the degree 
(or efficiency) of depletion (as well as the viability of different 
surface reactions) is explored in this study. Our in itial elemental 
abundances for Phase I are as in lBell etaD (120061) (see Table 1). 
We also ran some models where the initial abundances of both 
S and Si were either depleted or enhanced by a factor of 10 with 
respect to these values, reflecting the uncertainty of their degree 
of depletion onto dust. In Phase II, we follow the chemical evo- 
lution of the remnant core. For the hot core models, we simulate 
the effect of the presence of an infrared source in the center of the 
core or in its vicinity by subjecting the core to an increase in gas 
and dust temperature, up to T = 300 K. This increasing of tem- 
perature is based on the luminosity of the protostar by usin g the 
observational luminosity function of lMolinari et all (|2000). The 
models we use have been published befo re: for the model s de- 
scribing the hot cor e we refer the reader to lViti et all (l2004al) and 
iLerate et"aT1 ((2010) whil e for the models d e scribing the plateau 
we refer the r e ader to ILerate et aD d2008l) . ILerate et all d2010h 
and IViti et all (2004b). In both mo dels the presence o f a non 
dissociative C-shock (modeled as in Berg in et al.lll997l) can be 
simulated. If a shock is included in the model then sputtering 
also occurs and is faster than thermal evaporation. We have ran 
a total of 6 models for the hot core component and 2 models for 
the plateau component. For the hot core model we investigated 
the sensitivity of the chemical abundances to the degree of gas 
depleted on to the grains during the formation of the core; the 
branching ratios of surface reactions relevant to the formation 
of SiO and SiS; the initial abundances of Sulphur and Silicon; 
whether the gas is subjected to a non-dissociative shock during 
the hot core lifetime. For the plateau models, we only varied the 
initial Sulphur and Silicon abundances. 

We were able to reproduce the observed column density of 
SiO with most models. The only constrain we found was that 
the temperature of the gas must be at least ~ 100 K or, alterna- 
tively, must have undergone a shock. SiS, on the other hand, is 
difficult to produce: surface reactions (and subsequent evapora- 
tion or sputtering of the mantles) seem to be necessary. We find 
that the only models that succeed in reproducing the data are 
those where a percentage (even as small as 5%) of Sulphur on 
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the grains react with Si to form SiS: Figure [5] shows the column 
density of SiO and SiS as a function of time during Phase II of 
the hot core for models differing only in the mantle formation 
efficiency of SiS (i.e on how efficient Si bonds with Sulfur). A 
qualitative match with the observation can be achieved, at early 
times, by those models where only 5%-10% of Sulphur on the 
grains react with Si to form SiS, or at late times if a higher per- 
centage of Sulphur is involved in the formation of SiS. 

Note, however, that exact ages can not be derived from these 
considerations as the relationship of the time dependencies with 
the efficiency of SiS formation on grains will depend on the des- 
orption times of SiO and SiS. 

In conclusions, while it is possible to reproduce SiO in the 
gas phase (as well as on the grains), our models indicate that SiS 
is a product of surface reactions, most likely involving direct 
reactions of Sulphur with Silicon. 



6. Summary 

In Paper I we presented a line survey of Orion KL taken with the 
30m IRAM telescope. The sensitivity achieved allows to per- 
form a line confusion limited survey. Due to the wide frequency 
range covered and data quality we decided to present the line 
survey in a series of papers focused in different molecular fam- 
ilies. In this paper we presented the study of the emission from 
Silicon-bearing species, SiO and SiS, as well as their isotopo- 
logues and their vibrationally excited states. 

For the v=l state of SiO we have detected the J =2-1 line 
and, for the first time in this source, emission in the 7=4-3 tran- 
sition, both of them showing strong masering effect. The well 
known components of Orion (hot core, plateau, high velocity 
plateau, extended ridge, and compact ridge) contribute to the ob- 
served emission from SiO and its isotopologues whereas for SiS, 
29 SiS, and SiS v=l we have found emission from the hot core, 
the plateau, and a feature at 15.5 km s _1 . In order to check the 
origin of this feature we have observed the line profiles of 7=5-4 
SiO around IRc2 and of different molecules at an offset position 
(-15", 15") from IRc2. We conclude that this extended feature 
is seen around the position Aar=-7", A6-+7" (near the BN ob- 
ject) with a radius of ^5 arcseconds and probably arises from the 
interaction of the outflow with the ambient cloud. 

The physical parameters obtained for each Orion KL com- 
ponent are in agreement with those we can find in the Orion 
literature. For the feature at vlsr = 15.5 km s _1 we derive 
T/^200 K, n(H 2 )=5xl0 6 cm" 3 , and v /m(/ power ,„ ;e „.„- .=7.5 km 
s _1 . Column densities have been calculated with radiative trans- 
fer codes based on either the LVG or the LTE approximations 
taking into account the physical structure of the source and us- 
ing the new SiO-p-IL; and SiS-o-H2 collisional rates. The re- 
sults are provided as source averaged column densities. In this 
way, we obtain a total column density of (7.4+2.0)xl0 15 and 
(1.35±0.40)xl0 15 cirT 2 for SiO and SiS, respectively. 

We have derived several column density ratios which per- 
mit us to provide the following average isotopic abundances: 
28 Si/ 29 Si=26+10, 29 Si/ 30 Si= 1.7+0.6, 16 0/ 18 0>239, 18 0/ 17 0^2. 
We have also investigated the origin of the SiS and SiO emission 
in Orion KL by the use of a gas-grain chemical model and find 
that while SiO can be easily formed in the gas phase, SiS seems 
to be a product of grain surface reactions, most likely involving 
direct reactions of Sulphur with Silicon. 

The resulting vibrational temperature for SiS v=l in the fea- 
ture at 15.5 km s _1 is 500 K, larger than the kinetic temper- 
ature derived for this component indicating an IR pumping or 



a warmer component difficult to see in the lines of the ground 
vibrational state. 

Finally, we have derived upper limits for the column den- 
sity of non-detected molecules (Silicon-bearing species). For the 
detected species in other sources SiC, SiC2, c-SiC3, SiC4, SiN, 
SiCN, and SiNC the upper limits for their column density are 
1.3xl0 14 , 3.5xl0 13 , 1.3xl0 13 , 3.5xl0 12 , 6.1xl0 13 , 3.1xl0 13 , 
and 3.1xl0 13 cirT 2 , respectively. 
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Fig. 3. Observed lines (histogram spectra) and model (thin curves) of SiS. The dashed line shows a radial velocity at 15.5 km s~ 
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Fig. 4. Central panel shows a 20" x 20" map, centered at IRc2, of the 7=5-4 line of SiO; the right panel shows a comparison 
between the 7=5-4 line at an offset (-10",6"), and the SiS J=14-13 line towards IRc2, and the resulting model for each line (thin 
curves). The rest of the panels show the SiO J=5-4 line at different positions. 
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Fig. 5. Column densities of SiO (solid line) and SiS (dotted line) as a function of time from Phase II of the hot core model with 
an initial elemental abundance of lxl0~ 6 for Sulphur and 8x1 0~ 8 for Silicon. The different plots show different efficiencies for the 
formation of SiS on the grain, namely: 100% (top left); 0% (top right); 30% (bottom left); 10% (bottom right). 
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Appendix A: Density diagnostic 

In this section we study the influence of collisional rates on 
the modelling of the SiO and SiS lines observed in our line 
survey of Orion. Although most of the emission in Orion 
arises from regions of high volume density, the intensity of the 
high velocity wings apparent in the SiO and SiS line profiles 
arise from a region, the plateau, in which the density is not 
enough to thermalize the rotational levels of high dipole mo- 
ment molecules such as SiO (3.1 D) and SiS (1.73 D). Hence, 
collisional rates are needed for each species to derive the phys- 
ical conditions . For SiO previous calculations by Green and 
collaborators ( Bieni ek & Greenl 119831; iTurneret all 119921; see 
http://www.giss.nasa.gOv/data/mcrates#sio) have provided colli- 
sional rates for this molecule and temperatures 20-1500 K. These 
collisional rates have been re calculated using a new SiO-p-H2 
surface by iDavou & Balanca] d2006l) . Ho wever, no col l isional 
rate s were availab l e of S iS until recently: IVincent et al.l (2007) 
and lKlos&Liqud (f2008) have calculated the state to state colli- 
sional rates for the system SiS-He and SiS-H2, respectively. 

We have tested the influence of collisional rates on the mod- 
elling of the SiO and SiS lines observed in our line survey of 
Orion. Using the recent rates quoted above we have computed 
line intensity ratios between some transitions to provide tools 
to estimate the physical conditions in astrophysical sources. We 
have covered volumne densities between 10 2 -10 9 cirT 3 and ki- 
netic temperatures from 10 to 300 K including levels up to 7=20 
and 7=40 for SiO and SiS, respectively. In the following we will 
analyse the excitation conditions for 1— »0, 2— >1, 3— »2, 5— »4, 
12— »11, and 19— >18 for both molecules. 



A.1. SiO 

The line intensity ratios for SiO are shown in Figs. I A. 1 1 (ki- 
netic temperature between 10 and 100 K) and IA.2I (for a tem- 
perature of 300 K). In Fig. IA.1I we can notice that the inten- 
sity ratio Tb(3— >2)/Tfl(2— >1) is particularly sensitive to densi- 
ties of about 10 4 to 10 6 cm -3 , whereas for T B (5^4)/T B (l->0) 
it shows important variations in the density range 10 6 to 10 7 
cm -3 . Using the intensity ratio Ts(5— »4)/Tg(12— >11), we can 
explore density values around 10 7 -5xl0 8 cm 4 . For T/f=300 K, 
and T^(12— >1 1)/Tb(19— >18), we can trace densities above 10 5 
cm -3 , as seen in Fig. IA.2l For lower densities the 19-18 line will 
be very weak. 

To illustrate this point, we deduce from our observations the 
intensity ratio Tg(2— >1)/Tb(3— >2) of SiO, and find ^2.6 (see 
Table |2), which corresponds, depending on the column density 
and for T^=100 K, to densities between 10 5 5 and 10 7 cm 4 . Our 
results for the physical parameters of the different cloud compo- 
nents discussed in previous sections ar e given in Table [4) 

The collisional rates calculated bv lTurneretal]dl992l) have 
been used in the past to describe collisions of SiO molecules 
with H?. We have used through the paper the new rates of 
IDavou & Balanca! (2006). As a first step in modelling SiO emis- 
sion in warm clouds we have made a comparison between the 
results using both sets of collisional rates for a large range 
of physical conditions: T^=10-300 K, N(SiO)=10 12 -10 15 cm 4 , 
and n(H2)=l-10 9 cm 4 . Figure lA~3l sho ws the brightness temper- 
ature (T B ) ratio (predictions usi ng the Davou & Balanca, 120061 
rates over those obtained from iTurner et al.L 11992 rates) as a 
function of H2 density for different values of the SiO column 
density and temperature and for the six first transitions of this 
molecule. The plots in Fig. IA. 31 show that the difference in the 
predicted line intensities between both sets of collisional rates 



never exceed 40% (7=6-5 line), being always below 20% for 
all other rotational lines and kinetic temp eratures. The l owest 
temperature for the collisional rates of lTurner et al. ( 1992) is 20 
K; we have extrapolated these rates to obtain the corresponding 
ones at 1 K. For most transition s the predicted line intensities 
from the Dayou & Balanca ( 2006) rates are lower than those pre- 
dicted from ITurner et al.l (1 19921) rates for densities below ^10 5 
cm 4 . Although the differences in line intensities are not sig- 
nificant for the determination of the physical properties of the 
emitting gas in interstellar clouds, we have adopted i n our SiO 
calculations the new rates of lDavou & Balancal (|2006). 



A.2. SiS 

We have p erformed the same study for SiS using the rates cal- 
culated bv lVincent et al.l d2007l) . Figs. IA.4l and lA.5l show our re- 
sults. We observe the same trends as for SiO, but for lower den- 
sities. Thus, the Tg(3— >2)/Tb(2— >1) ratio increases in the den- 
sity range 10 3 to 10 5 cm 4 , the second one Tg(5— >4)/Tb(1 — >0) 
in 10 4 to 10 s cm 4 , and the last graphs of Fig. IA.4I show 
that T#(5— >4)/Tb(12— >1 1) vary mostly between 10 5 and a 
few 10 6 cm 4 . In Fig. [A3] we can see that for Tjc=300 K, 
T#(12— »11)/Tfl(19— »18) are useful for densities around 10 5 
cm 4 . Two of those lines have been detected in our survey, the 
transitions 5— >4 and 12— >11, with an intensity ratio of ^0.12. 
For T^=100K we can derive a volume density -10 5 5 cm 4 . 

Before 2007, no collisional data was available for SiS. 
Hence most authors have adopted the SiO collisional rates for 
Si S. As quoted above , the system SiS-He h as been studied 
bv IVincent et al.ld2007l) and lLique et al.l (|2008) have compared 
these rates with the correspondings for the SiS-H2 system. They 
have shown that, when scaled by the square root of the collision 
reduced mass, the SiS-He rates were a reasonable approximation 
to describe collisions with H2. We have updated the LVG code 
by adding the new SiS-He and SiS-H2 rates. We have calculated 
the line intensities and compared the resu lts to those obtained 
with rates for SiO from ITurner et al ] (1 19921) (multiplied bv a fac- 
tor 2 to take into account the larger geometrical size of the SiS 
molecule). Figure lA~6l shows the line i ntensity ratio (brigh tness 
temperature obtained with rates from Vincent et al., 2007 over 
those calculated with rates from ITurner et all 1 19921) as a func- 
tion of n(H2) for the six first rotational transitions. The discrep- 
ancies in the line intensities between the two cases do not exceed 
50%. As for SiO, we can notice that the differencies between the 
two sets of collisional rates are the same for all column densities: 
the line intensit y predicted with rat es for SiS are lower than those 
predicted from lTurner et all (11992) rates for densities up to ^10 3 . 
These differences are quite small and do not affect the determi- 
nation of the physic al parameters perfor med by authors that have 
used previously the ITurner et all Jl 992) rates in their models to 
interpret SiS observations. We tried also to describe SiS excita- 
tion b y collisions with the rate c oefficients for the system CS-He 
from Green & Chapman! (1 19781) with an appropriate scaling fac- 
tor. The predicted intensities show a reasonable agreement with 
those derived from SiS-He or SiS-H2 collisional rates. These re- 
sults show that with the actual calibration accuracy for the SiS 
observations, the determination of the properties of the emitting 
gas is not very sensitive to small differences in collisional rates. 
It is a rather curious result that the predicted line intensities with 
both sets of collisional rates, deduced from two different poten- 
tial energy surfaces, moreover, with different propensity rules, 
produce similar intensity ratios (when the scaling factor is ade- 
quately selected). In the determinatio n of the SiS co l umn d en- 
sities we have used the SiS-H2 rates (iKlos & Liauei 120081) al- 
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Fig. A.l. Intensity ratio between two transitions for SiO as a function of H2 density with the SiO column densities equal to 10 12 , 
10 13 , 10 14 and 10 15 cirT 2 and temperatures from 10 K to 100 K. 
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Fig. A.2. Same as|AJ]but for T = 300 K. 



though no significant differences are found when we used the 
SiS-He ones. 

In order to better quantify the effects of the use of scaled 
SiO-He rates in predicting SiS intensities, we have considered a 
case in which the lines are optically thin, for instance, n(H2)=10 5 
cm 4 , IV =40 K and N(SiS)=3.10 12 cm 4 . The predicted line in- 
tensities from SiS-He rates are T B (2->1, 4->3, 6^5)=0.08, 0.26, 
0.33 K. To obtain these results from scaled SiO-He rates we need 
a density of 7.1 4 cm 4 or to reduce the scaling factor from 2 to 
^1.5, i.e., an error lower than 2 in density. 

We have compared predictions from the SiS-He rates with 
those obtained from the SiS-H2 rates and we have found not sig- 
nificant variation between both. 
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Fig. A.3. Comparison of the line intensities predicted using two different sets of collisional rates (see text) for different lines, column 
densities and densities. 
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Fig. A.4. Intensity ratio between two transitions for SiS as a function of H2 density with the SiS column density equal to 10 12 , 10 13 , 
10 14 and 10 15 cm" 2 . 
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Fig. A.5. Same as|A3]but for T = 300 K. 
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Fig. A.6. Same as|A3]for SiS with the SiS column density equal to 10 12 , 10 13 , 10 14 and 10 15 cm - 
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Fig. B.l. Integrated intensity of the 7=5-4 transition of SiO at different velocity ranges (indicated at the top of each panel). In some 
panels, the integrated intensity of the maps has been multiplied by a scale factor (indicated in the panels) to maintain the same color 
dynamics for all maps. The step in integrated intensity (J T^dv) is 20 K km s" 1 from 25 to 305 K km s _1 and the minimum contours 
are 5 and 15 K km s _1 . See text (Sect. [2]) for main beam antenna temperature correction. 
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Fig. B.2. Observed lines (histogram spectra) and model (thin curves) of SiS and SiS v=l. The dashed line shows a radial velocity 
at 13.5 km s _1 . 
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Fig. B.3. Observed lines of different molecules at a position (-15", 15") offset IRc2. The dashed line shows a radial velocity at 15.5 
km s _1 . 
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Fig. B.4. Observed lines (histogram spectra) and model (thin curves) of SiN. 
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Table C.l. SiO and SiO isotopologues velocity components. 



Species 




Wide component 




Narrow component 






\ LSR (km s" 


') Av(kms-») 


Tmb (K) 


v iS /i (km s _1 ) 


Av (km s 1 ) 


Tmb (K) 


SiO 2-1 


9.4±0.5 


44±2 


5.19 


8.45±0.13 


19.1+0.5 


13.3 


SiO 3-2 


10.2±0.3 


40.1+0.9 


15.5 


8.4±0.2 


14.9±0.7 


10.8 


SiO 4-3 


12.0+0.3 


41.7+0.9 


26.6 


8.5±0.2 


15.8±0.8 


18.8 


SiO 5-4 


11.91+0.13 


40.4±0.4 


27.1 


7.30+0.11 


14.8±0.4 


18.4 


SiO 6-5 


9.7±0.3 


40.8+0.5 


46.1 


7.7±0.3 


13.3±0.8 


15.6 


M SiO 4-3 


7.9+1.1 


46+6 


2.88 


9.5+0.3 


15.5+1.3 


5.85 


30 SiO 4-3 








8.8±0.2 


19.7±0.5 


5.15 


30 SiO 5-4 








10.23±0.15 


21.0±0.6 


6.19 


30 SiO 6-5 


7.8±0.5 


46+2 


4.35 


8.6+0.2 


16.8±0.7 


8.50 



Notes. Velocity components of selected SiO and SiO isotopologues lines derived from two Gaussian fits. 
Table C.2. SiS velocity components. 



Species 




Wide component 




Narrow component 






Vlsr (km s' 


- 1 ) Av(kms-') 


Tmb (K) 


v LSR (km s~') Av (km s" 1 ) 


Tms (K) 


SiS 6-5 


8.8±0.5 


24.0+0.9 


0.39 


15.4±0.3 7.0+0.5 


0.23 


SiS 8-7 


9.6±0.4 


24.9±1.0 


1.03 


15.7±0.5 6.3±1.0 


0.49 


SiS 9-8 


8.7±0.4 


24.6±0.7 


1.33 


15.7±0.2 6.6+0.5 


1.15 


SiS 15-14 








15.7+0.3 12.3±0.9 


2.69 



Notes. Velocity components of selected SiS lines derived from Gaussian fits. 
Table C.3. Molecular abundances 



Region 


Species 


X (This work) 


X' 


X s 


X y 






(xl0~ 8 ) 


(xlO- 8 ) 


(xlO- 8 ) 


(xl0~ 8 ) 


Extended 


SiO 


0.05 


<0.05 






ridge 1 


SiS 










Compact 


SiO 


0.20 


0.80 






ridge 2 


SiS 










Plateau 3 


SiO 


2.20 


0.80 


1.10 


23.5 




SiS 


0.30 






0.40 


High velocity 


SiO 


2.20 




4.50 




plateau 4 


SiS 










Hot 


SiO 


0.20 


0.60 






core 5 


SiS 


0.30 








15.5 km s-' 


SiO 


1.00 








component 6 


SiS 


0.70 









Notes. Derived molecular abundances and comparation with other works. 

(1) Assuming AT(H 2 )=7.5xl0 22 cnr 2 . (2) Assuming A'(H 2 )=7.5xl0 22 cnT 2 . <3) Assuming A r ( H2)=2.1xl0 2 3 cm' 2 . <4) Assu ming N(H 2 )=6 2xlCP 
cmr 2 . (5) Assuming jV (H 2 )=4.2xl0 23 cm" 2 . (6) Assuming A'(H 2 )=1.0xl0 23 cm 4 . <7) From ISutton et all {T995). <8) From iPersson et all j2007h . 
(9) From lZiurv1 lfr988). 



